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Abstract 

We investigate the pseudo-Casimir force acting between two charged 
surfaces confining a single polyelectrolyte chain with opposite charge. 
We expand the exact free energy to the second order in the local elec- 
trostatic field as well as the replicated polymer density field around 
the mean-field (saddle-point) solution. The quadratic terms lead to 
a fluctuation interaction that is partly due to the (thermal) Casimir 
effect for the confined electrostatic field and partly due to the pseudo- 
Casimir effect due to the confined replicated polymer density field. 
We study the intersurface separation dependence of both effects and 
show that the pseudo-Casimir effect leads to a long range attraction 
between the surfaces that decays with an anomalous algebraic expo- 
nent of ~ 1.7, smaller than the standard exponent of 2 in the case of 
Casimir interactions. 

1 Introduction 

Charged polymers have received much attention over the years for their 
indisputable importance in industrial applications as well as for their role 
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in fundamental colloid science (for comprehensive introductions see ^|). 
Apart from the study of polymers in the bulk, their interactions with sur- 
faces, as well as polymer-mediated interactions between surfaces have been 
studied on different levels [||, |J. Though polyelectrolyte interaction with 
charged surfaces is quite similar to the general interaction of polymers with 
neutral surfaces ||, the long range nature of the Coulomb potential in- 
troduces additional features and difficulties that make these problems even 
more difficult to handle. 

Most of the work on inhomogeneous polyelectrolytes was done on the 
mean-field level and the effect of thermal fluctuations has not been consid- 
ered at all. The mean-field picture of the behavior of a polyelectrolyte chain 
confined between two oppositely charged surfaces that emerged from previ- 
ous studies [0, ^] (see Fig. [I]) lets one think that at small intersurface sepa- 
rations the polyelectrolyte mediated interations are quite similar to those in 
the case of neutral polymers and decay algebraically with separation || ||] . 
At a critical value of the intersurface spacing there is a transition of the 
confined polyelectrolyte chain from a monomodal to a bimodal monomer 
density distribution. This means that at sufficiently large intersurface sepa- 
rations the polymer chain is on the average partly electrostatically adsorbed 
to both surfaces and the remaining part of the chain bridges the region be- 
tween the surfaces. The consequence of these bridging interactions is an 
exponentially decaying attraction between the surfaces that persists for all 
separations greater than the critical value 0, ^. In this respect the charged 
chain between oppositely charged surfaces behaves on the mean-field level a 
lot like a neutral chain with a finite adsorption energy to the two bounding 
surfaces. 

Recent theoretical studies || however point to conclusion that thermal 
fluctuation effects in Coulombic systems, ignored on the mean-field level, 
can hardly ever be considered as negligible. In the limit of small electro- 
static coupling it has been known for a while that thermal fluctuations in an 
inhomogeneous Coulomb fluid lead to zero order van der Waals (or equiva- 
lently Casimir) forces |TT| , ^]. More recent studies of thermal effects in soft 
matter in general Jl^, are very much consistent with these findings. It 
has been established that thermally driven fluctuations of order parameters 
of confined systems also lead to long range pseudo-Casimir forces |10], quite 
similar in nature to ordinary Casimir forces except that instead of being 
due to electromagnetic field fluctuations they are a consequence of a general 
order parameter fluctuations (thus pseudo-Casimir interactions). In view of 
this we suspect that similar phenomena should also show up in the study of 
fluctuation effects of confined inhomogeneous polyelectrolytes. 

In order to prove this conjecture we write the partition function for a 
polyelectrolyte chain confined between two oppositely charged planar sur- 
faces in such a way that we were able to expand it to quadratic order in 
the relevant order parameters, in this case the local electric field and the 
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replicated polymer density field (see below), around the mean- field profile 
and evaluate on this level the effect of thermal fluctuations. We note that 
thermal fluctuations decouple into two separate contributions: the Casimir 
contribution due to electric field fluctuations and the pseudo-Casimir con- 
tribution due to the replicated polymer density field fluctuations. We show 
that the effect of thermal fluctuations depends on only one coupling param- 
eter olmf = jf, where a is the surface charge density, r is the line 
charge density of the polymer chain and I is the Kuhn's segment length. 
Depending on the value of this parameter the behavior of the system can 
be successfully described on the mean-field level [olmf 3> 1) or has to be 
ammended by taking into account the fluctuation terms [olmf ~ !)■ The 
most important consequence of the fluctuation effects seems to be the long 
range pseudo-Casimir force, with an anomalous algebraic dependence on the 
intersurface separations with an exponent 1.7, due to the confined fluctua- 
tions of the polymer density field that decays slower than the usual Casimir 
(or equivalently van der Waals) force. 

The plan of the paper is as follows. In order to investigate the fluctuation 
effects we first write the free energy of a confined charged polymer chain 
between two surfaces in the form of an 0{n) scalar field theory. From 
here the free energy is obtained via the standard n — > replica trick 14]. 
Since the mean-field is given by the saddle-point, we then expand this field 
theory to the second order around the saddle-point in both the electrostatic 
field and the replicated polymer density field. Calculation of the quadratic 
corrections to the mean-field free energy that can be performed exactly 
gives the Casimir and pseudo-Casimir interactions mediated by the confined 
electrostatic and confined replicated polymer density fields. We investigate 
the dependence of the regularized fluctuation interaction free energy on the 
intersurface separation in the limit of small and large spacings compared to 
the critical separation. 



2 Mean-field theory 

As a point of departure we shall take a model confined polyelectrolyte sys- 
tem, composed of two surfaces with a specified surface charge density (a) 
and a confined, oppositely charged single polyelectrolyte chain of length N£, 
described with an Edwards hamiltonian 

H[r(s)} = | / r 2 ( S ) ds + y M J M V(r(s),r(s')) dsds' . (1) 
J o 

We assume that the intersegment interactions are given by the unscreened 
Coulombic form V(r, r') = r 2 /47reeo |r — r'|, where r is the charge per unit 
length. This is a bit artificial since we completely ignore the presence of 
other mobile ionic species such as counterions and salt ions, but represents 
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a clean, tractable model system where fluctuation effects can be studied in 
detail. 

Performing now the Hubbard-Stratonovich transformation for the pair 
interaction potential, going to the grand canonical ensemble with a fixed 
chemical potential [/, for the monomers and adequatly removing the closed 
loop polymer configurations from the partition function via the n — ► limit 
of an 0{n) scalar field theory with the field ip = (ipi, ip2 ■ ■ ■ "0n)> we remain 
with the de Gennes - des Cloizeaux representation for the free energy of a 
charged polyelectrolyte chain Jll 



F(jji,P) =Jbm o ilogZ(/x J /3)=lim n ^ ilog/^(r)PV'(r) exp (-/? / d 3 r H^r), ^(r)]) 

(2) 

where 

/ 9H[0(r),^(r)] = if Ei^fi^ £, + f/^o W) 2 + z/3r E^Vf 

(3) 

In writing Eq. [2] we omitted the determinant det (/3eeoV 2 ) from the de- 
nominator. As is usual on the one-loop expansion level, this term is exactly 
cancelled by the zero order van der Waals term that one should add to the 
final total interaction free energy. This point is thoroughly explained in 



[16, 11]. See also below. 

The path from here will be to obtain the mean-field solution of the model 
defined by Eq. ^ via the saddle-point of the 0(n) scalar field theory along the 
lines of |l8| and then to evaluate the contribution of quadratic fluctuations 
around the saddle-point to the thermodynamic properties of the system. 

Evaluation of the saddle-point from the field hamiltonian Eq. || is triv- 
ial. One notices immediatley that all the tpi satisfy the same saddle-point 
equation and thus ipi — * ip. With this one can write the dependence of the 
hamiltonian Eq. || on n explicitly [19] and the limit n — ► can be obtained 



straightforwardly. One remains with 

(3Fo = - Hm i log Z = i / d 3 r ( f ( ViP) 2 +/i^ 2 + f3ee ( V0) 2 + ifir <f>^ 

(4) 

where we introduced the mean polymer density field tp so that the local 
monomer density is p = ip 2 and the local electrostatic mean potential tp = i<p. 
In the planparallel geometry considered here all the fields depend only on 
the transverse coordinate z, and the inhomogeneities are confined to this 
direction only ip = ip(z), (p = <p(z). The saddle-point equations of the 
hamiltonian Eq. [3| can thus be written in the form 

I 2 d 2 ip(z) ,. , nn , . . d 2 ip(z) 2 , 

g~ ffi- = + 2pT<p{z)i>(z) , -ee ^^=T^ 2 {z). (5) 
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These mean-field equations are exactly equivalent to those derived previ- 
ously (7|, ||] if one takes into account that the polyelectrolyte chain is the 
only mobile charge in the system (no counterions and no salt). The first 
equation of Eq. || is the equation for the density field ip(z) of the polymer 
in an inhomogeneous external field <f(z) while the second one is a Poisson - 
Boltzmann equation for the mean electrostatic potential f(z) of a charged 
polymer chain with charge density p(z) = Tip 2 (z). 

The second equation of Eq. [| can be solved explicitly and its solution can 
be manipulated JjJ to yield a limiting form valid close to the boundaries at 
z ~ ±a where the potential is largest 

\z - z'\^ 2 {z') dz' ~ -t\z\ + 0{z 2 ). (6) 

-a 

This expression can be obtained rather straightforwardly by Taylor expand- 
ing the solution around z = ±a. In the opposite limit, z ~ 0, we can derive 

/+a 
\u\ip 2 (u)du - Tip 2 (0)z 2 + C(z 4 ). 
-a 

Since the mean-field is largest close to the boundaries we will use the approx- 
imate expression valid strictly only close to z = ±a in the whole interval. 
This approximation works extremely well everywhere |7j except close to the 
origin where z is small anyway. The reason for introducing this additional 
but inessential approximation is that it also helps in transforming the first 
of Eqs. H from a non-linear Landau-Ginzburg type equation into a linear 
equation, which can be solved analytically and explicitly in terms of Airy's 
functions 0]. 

The electroneutrality demands that at the bounding surfaces (z = ±a) we 
have — eeo§^7 = r 3r = 0*j where n is the boundary surface normal. Because 
of the impenetrability of the boundaries to the polymer chain we should 
also have ip{z = ±a) = 0. The unnormalized solution of the first mean-field 
equation Eq. [5| can now be obtained explicitly in the form 

ip(x) ~ Ai(y - x)Bi'(y ) - Bi(y - x)Ai'(y ), (7) 

where we introduced the dimensionless variables x = X^ 3 z, xq = A^ 3 a, 
yo = 6//A B 2//3 /£ 2 , A# = 7^2- and Ai(x) and Bi(x) are the standard Airy 
functions. The dependence yo = yo(xo) is obtained from the vanishing 
density field boundary conditions at the two bounding surfaces. We derive 
0] that yo(xo) grows linearly with xq for large xq, and therefore u$ = yo — ^o 
asymptotically approaches a constant Uq° ~ —2.34, while for small xq, yo is 
negative and behaves approximately like yo rs — (it/2xq) 2 . 

The main consequence of the mean-field equations is that for very small 
intersurface separations the chain is desorbed with a monomodal density 
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(C) 




Figure 1: Schematic representation of the mean- field solution. For sepa- 
rations between the surfaces smaller then the critical, the monomer den- 
sity distribution is monomodal (A). For separations larger then critical the 
monomer density distribution is bimodal (B), corresponding to partial ad- 
sorption of the chain to the oppositely charged surfaces. Parts of the chain 
that are not adsorbed act as bridges between the surfaces (C). 



distribution @, §]. In this regime of intersurface separations the interac- 
tions between the surfaces are repulsive, decaying albebraically with plate 
separation. At a critical value of the intersurface separation xo = 1.986 
the electrostatic attraction causes the chain to adsorb to both surfaces cre- 
ating a polymer bridge leading in its turn to bridging attraction between 
the surfaces. This attraction, though exponentially small, persists for all 
separations between the surfaces larger than the critical. The intersurface 
pressure is thus repulsive for xo < 1.986 and attractive for xq > 1.986. For 
details see |7]]. 



3 Fluctuations 

We now turn to the contribution of fluctuations around the saddle point 
to the free energy. First we expand the hamiltonian Eq. |3| to the second 
order in deviations from the mean-fields Eq. |5[ Our treatment of the fluc- 
tuation effect will thus be based on a quadratic (one-loop) expansion. The 
second functional derivatives of the field hamiltonian can be assembled into 
a Hessian of the fields of Eq. |3| 

_ [ -fv 2 + n + 2i(3T<t) 2i/3r^„ 
[ 2i(5r^ n -pee V 2 

where cp and tp n are the solutions of the mean-field equations Eq. [|. While 
integrating out the mean replicated polymer density field and the mean 



(8) 
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electrostatic field fluctuations we notice that the functional integral over ipi 
fluctuations can be evaluated directly, giving an explicit dependence of the 
fluctuational partition function on n. This is similar to the general case 



treated by Emery [19]. The limit n — > can now be evaluated explicitly 
leading to the following contribution of the fluctuations around the mean- 
fields in the quadratic order to the free energy 

0F 2 = lim i log Z 2 = ^lndet (- V 2 ) + ±lndet (- V 2 + £/z + 12ffa(r)) + 
+ 12^ J J d 3 r dV g(r,r') £ (r,r') ^(r)^(r'), (9) 

where <£>(r) and ijj(r) are the solutions of the mean- field equations Eq. ||. 
The two Green's functions Q(r,r') and £/o( r > r ') are the functional inverses 
defined as 

(-V 2 + |M+12^(r))g(r,r') = o 3 (r-r') 

(-V 2 )g (r,r') = o 3 (r-r'). (10) 

The first term in Eq. ^ comes from the confined electric field fluctuations 
and is in fact the standard zero-order Casimir interaction energy |2(| (see 
below). The second and the third terms are due to the confined polymer 
density field fluctuations and thus correspond to the pseudo-Casimir inter- 
actions. The free energy Eq. ^ is obtained already after the cancellation 
referred to in the discussion following Eq. || and thus represents the total 
interaction free energy including the zero order van der Waals contribution. 
The latter is evaluated for a slab of thickness 2a with boundaries impene- 
trable to electrostatic fields. 

Let us first deal with the two functional determinants in Eq. ^. Since the 
mean-field solution depends only on the transverse coordinate z we need to 
investigate the functional determinant of an operator of the type — -py + 
V(z) that can be evaluated in two different ways. One can first of all find 
the eigenvalues A n of this operator with the boundary condition f n (z = 
±a) = 0. The functional determinant comes out as the product of these 
eigenvalues U n \ n . The other approach would take into account the van 
Vleck identity [21] involving fo(z), the eigenfunction with zero eigenvalue 
and with boundary conditions fo(z = —a) = and f${z = —a) = 1. The 
logarithm of the functional determinant can thus be written equivalently 
(up to an irrelevant - in this context - additive constant) in two different 
ways |2l| 

In det (-^ + V{z) \ = InAn = In [f (z = a)} . (11) 

In what follows we will use both representations together with the Fourier 
decomposition in the (x, y) plane and will find both equally inadequate 
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when approaching the critical intersurface separation xq = 1.986, either 
from above or below, where the mean polymer density field goes through a 
monomodal to bimodal transition. 

Of the first two terms in Eq. ^ the first one is trivial to evaluate either 
via the van Vleck relation or the eigenvalue product representation. First 
of all we introduce the Fourier decomposition in the (x, y) plane with a 
wave vector Q, since we only need to take into account the boundaries in 
the z direction. This leads to the standard zero-order (thermal) Casimir 
interaction in the form 

\k B T E Q In det = Skf. f°° QdQ In ( coshQa sinhQa \ ^ {u) 



where S is the total area of the interacting boundary surfaces. 

The second term is trickier. First of all we again introduce the Fourier 
decomposition in the (x, y) plane with a wave vector Q and invoke the 
approximate solution Eq. |(| Using again the van Vleck identity we derive 

\k B T E Q In det (-£ + Q 2 + £/x + 12§^(z)) = ^f- j °° QdQ In (C(Q, a)S(Q, a)) , 

(13) 

where we decomposed the determinant of the operator into two separate 
contributions stemming from symmetric and antisymmetric modes 

C(Q, z) = Bi'(y(Q))Ai(y(Q) - X B /3 z) - Ai'(y(Q))Bi(y(Q) - xf 'z) 
S(Q, z) = Bi{y{Q))Ai{y{Q) - X^z) - Ai(y{Q))Bi{y{Q) - A^)(14) 

where y(Q) = X B 2 ^ (q 2 + ^ ) while Ai{x) and Bi(x) are the standard Airy 



functions |g2[ . The Airy function solutions are of course a consequence of the 
approximation Eq. |B] for the mean electrostatic potential. The free energy 



Eq. 13 still needs to be properly regularized by subtracting the infinite bulk 
and surface terms, see below. 

In what follows we introduced also these dimensionless variables: q 2 = 
X B 2/3 Q 2 and u(q,x,x ) = y(q,x )-x = q 2 +y (x )-x = u (x ) + q 2 + x -x. 
Here yo(xo) is the largest zero of the symmetric function C(yo(xo),xo) which 
corresponds exactly to the zero of the mean-field solution Eq. [?] at the two 
bounding surfaces M . 

The third term in Eq. || appears to be the most difficult one to evaluate 
but fortunately it is possible to prove that within the approximation Eq. || its 
contribution to the fluctuation free energy Ti is in fact negligible compared 



to the other two, i.e. Eqs. 13 and 14. Taking into account the approximate 



form for the mean electrostatic potential Eq. y the last term in Eq. |9j can 
be evaluated explicitly leading to 

12 ^F / / ^ dV G{r,r ' ] G ° ir,r ' ] ^ (r ^ (r/) = 7^) f +a dAAG^z',Q = 0), 

(15) 
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where we took into account that the mean-field solution ip(r) depends only 
on the transverse coordinate z and that again we can introduce the Fourier 
decomposition in the (x,y) plane with a wave vector Q. The magnitude 
of the coupling contribution Eq. turns out to be much smaller than 



the other two, Eqs. [L^, 13, and can thus be safely ignored in numerical 
computations. 

The fluctuation free energy, T 2 which is thus composed only of the contri- 
butions from the fluctuational determinants of the electrostatic field and the 
replicated polymer density field, still contains the divergent bulk and surface 



parts together with the interaction terms [11, ^C|. T 2 is thus formally infi- 
nite. This infinity can be regularized by subtracting the infinite bulk and 
surface contributions |2(J and thus in forming the regularized fluctuation 
free energy 

m(a) = F 2 (a)-F?{a). 

Here J-2°(a) stands for the form of the fluctuation free energy for large values 
of the argument but evaluated at a finite a |20| ], 

Using now the standard form for the mean-field pressure |8|] and adding to 
it the fluctuation contribution, one can write for the total pressure between 
the boundary surfaces 



p(a) = pmf + Pfl [5<t>] + Pfl [ 

= ik B Tptf(z = 0)-^$$+p FL (a), (16) 

where the first term on the r.h.s. derives from the mean-field solution, the 
second one derives from the Casimir term (stemming from the electrostatic 
field fluctuations, first term in Eq. |9] and the third one derives from the 
pseudo-Casimir terms (stemming from the replicated density field fluctua- 
tions, second term in Eq. with 

1 ( dW 2 {a) \ 

PFL{a) = -s{-War)- { ] 

The next section deals with the evaluation and regularization of the fluctu- 



ation terms in Eq. 16 



4 Results - large intersurface separation 

We will first have to regularize the divergent expression Eq. [l^ by subtract- 
ing the divergent terms explicitly in the limit of large xq (to be quantified 
below). When xq is large enough the argument y = u + xq of the Airy func- 
tions in Eq. [l4| becomes large too and we can use the asymptotic expansion 
formulas |^2| for the Airy functions to explicitly find the divergent terms in 
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Figure 2: The regular part of the fluctuational free energy F = W2 in 
the limit of large separations (see main text) plotted with respect to the 
dimensionless separation xq. The region around the critical separation xq = 
1.986, where the approximations are doubtful, is marked with a box. On 
the inset there is the same plot in the logarithmic scale (log | w I) • One can 
read from the inset that the dependence of the free energy on the separation 
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is algebraic with the exponent around 1,7 (W2 oc x ' ). The constant F* 
is defined as: F* = — 167r B . 

•^ 7 2°( a )- The regularized free energy is then 

mixo) = ^i r du 

47r Ju 

j / (Si' (it + x )Ai(u) - Ai'(u + x )Bi(u)) (Bi(u + x )Ai(u) - Ai(u + xq)B%(u)) 
I (Bi' L (u + x )Ai(u) - Ai' L (u + x )Bi(u)) (Bi L (u + x )Ai(u) - Ai L (u + x )Bi(u)) 



where Ail, Bzl, Ai' L and Bi' L are the lowest order asymptotic expressions 
for the Airy functions |22j| . The next order in the asymptotic expansion of 
the Airy function is AiL2{t) ~ — 0.07(|t 3 / 2 ) _1 Ai^(^) and in order for it to 
be negligible, t has to be large enough, typically larger then 1.5. This leads 
to the condition uo(xq) + xq > 1.5, which yields roughly that xo has to be 
larger than ss 3 if our calculation is to be accurate to within a few percents. 
For smaller values of the intersurface separation the form of the regularized 
interaction free energy Eq. [l^ does not work. 

Fig. H shows the dependence of the dimensionless fluctuational free energy 
W2 on the dimensionless separation xq. As already stated the regularization 
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procedure Eq. Il8] can only be used for sufficiently large intersurface spac- 



ing. The log-log plot (see the inset of Fig. g) reveals that asymptotically 
the fluctuational contribution with an algebraic decay prevails. From the 
inset of Fig. || we read that the free energy decays asymptotically as Xq 1,7 . 
The algebraic decay of the free energy and concommitantly of the pressure, 
indicates that the fluctuation induced force has a much longer range than 
the mean field force which decays approximately exponentially ffl, |8| . Also 
it has a longer range than the zero-order Casimir (or equivalently the zero 
order van der Waals) interaction free energy that decays 

There is also no obvious connection between this result and the behavior 
of the fluctuation part of the interaction free energy of a system composed of 
unconnected counterions treated on the same level of approximmation |23] . 



In that case the fluctuation free energy decays as x 2 Iiixq, which is much 
faster than x^ 1 ' 7 . 



5 Results - small intersurface separation 



For small plate separations the above procedure cannot work because in 
that case we can not write down explicitly the asymptotic form of the Airy 
functions that would be valid in the whole range of integration in Eq. [Dj| . 
We found out that the best way to approach the evaluation of the fluctuation 
determinant in this limit is via a perturbation analysis using the eigenvalue 



representation of the functional determinant Eq. 11 



The perturbative treatment of the evaluation of Eq. [D| can be introduced 
as follows. The differential equation we have to solve is 



d 2 

-j^2+<l 2 + yo(^o) + \x 



x \ )ip = . 



(19) 



The last term can be treated as a linear perturbation potential V(x) = 
\x — xq\. What we will do now is to evaluate the determinant of the oper- 



ator in Eq. 19 in the form of a perturbation expansion in V(x) up to and 
including the first order. The order zero of the perturbation expansion of 



the determinant of the operator Eq. IS gives 



det 



dx 2 



+ q +yo 



sinh2 v / g 2 + y x 

iVv 2 + Vo 



(20) 



The final integration over the wavevector Q has to be divided into two parts 
because yo is negative, leading to 



W 2 ° 



Sk B T\^ 3 

167TXQ 



sinw C(3) 
wdw m — 

ij w 4 



(21) 
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Figure 3: The regular part of the fluctuational free energy F = W2 in 
the limit of small separations. The region around the critical separation 
xq = 1.986, where the approximations are doubtful, is again marked with a 
box. The zero - order approximation is shown dashed and the result up to 
the first order solid. The error bars represent the maximal possible error of 
the first order result. The inset represents the log- log plot of the free energy, 
wherefrom we deduce that for small separations the fluctuation free energy 
scales . The constant F* is defined in Fig. |2|. 



To obtain the first order term in the perturbation expansion we shall 
use the product of the eigenvalues representation of the determinant of the 
operator which leads to the following first order form of the free energy 

,2/3 



T\ = 3 / qdq £ InAjv (22 

^ Jo N 



where X\r represents the first order correction to the eigenvalue of the dif- 



ferential equation Eq. IS. This correction can be evaluated by standard 



formulas of the perturbation theory p4| ] leading to the following result 

C 2x I I ■ 2 ( Nit \ j 

a _. <N\V(x)\N> _. Jo \x-x Q \sm [^x) dx 

Kn ~ Xn <^v|iv> - x * r^^f^H, " N ' 2 

(23) 



! _ , <N\V(x)\N> o Jq F-xo| S m ^xjax _ 



so the sum of the first order corrections to the log of the determinant of the 
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operator Eq. 22 is 



N 



N 



n N 



N 2 



i+t^fw+f; 



(24) 



Since the first term in this expression has already been calculated, we now 
focus on the second one. This term is obviously divergent. If we calcu- 
lated the bulk contribution, there would be an integral over N instead of 
the discrete sum, which is a consequence of the boundary conditions. The 



regularized free energy is thus [2C] obtained as the difference between the 
sum and the integral of the same expression. After some algebra the result 
of this manipulation is 



! = Sk B TX 2 ^x 
2 16vr 



+ 



dw ^7T (coth TTW — COt TTW — 1 
-3 dw ( 7T (coth TTW 



,1 + 10 „ 1 N 

In (- iarctan — 

1 — 10 w / 



+ 



1) + 2 ( arctan— 

y w 



5A]v < 



1 
AA 



< N\V(x) 2 \N > - (< N\V(x)\N >) z 



(26) 



Even without really calculating the second order we can nevertheless say 
something about its magnitude. The maximal difference between the first 
and the second order perturbation terms can be obtained in standard nota- 
tion as pl 



where AA is the minimal difference between A^v and any other Ajvf- A short 
calculation then yields %- < j^g. As expected, the reliability of the first 
order result drops as we move towards larger xq i.e. towards the mean field 
phase transition from bellow. 

Fig. ^ shows the first and the second perturbation order of the regularized 
fluctuation free energy at small intersurface separations. The pressure ob- 
viously changes sign for sufficiently large separation (xo larger than ~ 1.68 
on the Fig. ||). For small enough separations the free energy scales ap- 
proximately as Xq 2 which is the same scaling form as for the mean-field 
part. It seems to us that the change in sign of the pressure close to the 
transition point xq = 1.986 is not due entirely to artifacts of the first order 
perturbation theory for sufficiently large intersurface separations, but could 
concievably survive all higher order perturbation terms. Unfortunately all 
approximation schemes that we could think of broke down close to the tran- 
sition point and we thus have to leave the important question about the 
behavior of the fluctuation component of the pressure close to the transition 
unanswered. 



(25) 
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Figure 4: The mean field pressure (dot-dashed line) and the total pressure 
{aMF = 4, thick line and aMF = 100, thin line) - between two plates as 
a function of the separation in the logarithmic scale. The uncertain region 
around xo = 1.986, where the approximations are doubtful, is marked by a 
box. The box contains only the mean-field results which can be evaluated 
for the whole range of xq values. One can clearly observe the long range tail 
stemming from the thermal fluctuations of the monomer density field. Inset 
shows the total pressure - mean field plus fluctuational - between two plates 
as a function of the separation for different values of parameter aMF = 
4, 100. A large deviation from the mean field result is obtained for aMF = 4 
whereas fluctuations hardly matter for aMF = 100. The constant p* is 

2 

defined as: p* = 

6 Discussion 

We have described all the contributions to the fluctuational free energy of 
the system. To evaluate the fluctuational part of the pressure between the 
two surfaces the derivatives of the regularized free energy with respect to xq 
have been performed numerically. The total prerssure can thus be obtained 
by simply taking the sum of mean-field and fluctuational contributions Eq. 
|l6[ Fig. [| shows the dimensionless total pressure as a function of the 
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dimensionless separation xq which can be written, Eq. 16, as 
p(x ) = pmf +Pfl[S4>] +Pfl[ 



k B T\ B £ 2 a J2 k B T({3)\ B .1/3 d fW 2 (x )\ 

— — yo (z o) ^ (0) - 64 ^3 " A B 



P 



( \ 72 /n\ C(3 /(go) 

yo{x )ip (0) - 3 + 

4 OiMF xfi a M F 



(27) 



where /(xq) is a function of xo only , V>(0) is the normalized dimensionless 
^ QMF = j7and/= fcBT i 6 B " MF = The parameter qcmf depends 
on the charge density on the plates cr, the linear charge density on the 
polyelectrolyte chain r and on the segment length of the polyelectrolyte I. 
If we put in experimentally reasonable values of a ~ eo/1, 35nm 2 then omf 
is about 100 for DNA (a relatively stiff polyelectrolyte) and around 4 for 
hyaluronic acid (a relatively flexible polyelectrolyte). Obviously the larger 
aMF the less important will the fluctuation effects be and the more accurate 
will be the mean-field result. 

Fig. H] shows the complete pressure Eq. |27] as a function of the dimen- 
sionless separation for different values of the coupling parameter aMF- As 
already stated we were unable to analyse the behavior of the interaction 
pressure close to the mean-field transition point. Our failure could indicate 
that close to the transition point our harmonic approach to the fluctuations 
breaks down. This would not be totally unexpected since we are dealing 
with a continuous transition in the distribution of the polyelectrolyte den- 
sity that apparently has all the attributes of a second order phase transition, 
with large fluctuations close to the transition point. We will however not 
pursue this line of thought here. 

It is clear from Fig. || that the effect of the fluctuations on the intersurface 
force can be profound. First of all the fluctuations modify the force at small 
separations making it less repulsive if compared to the mean-field case. Since 
both the mean-field as well as the fluctuation contribution to the free energy 
scale as Xq 2 (but with a different sign) it appears that the effect of the 
fluctuations is simply to renormalize the magnitude of the pressure but not 
its scaling form. If one would thus try to fit the complete p(a) with pmf{o) 
one would have to introduce a smaller effective charge on the surfaces for 
the fit to make sense. 

At large separations the behavior of p(a) is completely different from 
the behavior of pmf{o)- What the fluctuations do is that they provide 
a long range algebraic tail to the interactions, which overwhelms the ex- 
ponentially screened mean-field attraction at large intersurface separation. 
The algebraic fluctuation tail decays with an anomalous exponent of 1.7, 
that makes the polyelectrolyte mediated fluctuation interactions even longer 
ranged than the standard van der Waals interactions. This in itself is one of 
the main conclusions of our work. It is also this long ranged tail that is prob- 
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ably most amenable to experimental observation and verification. Since the 
fluctuation effects are stronger for smaller olmf we predict that they should 
be easier to measure for very flexible polyelectrolytes, such as hyaluronic 
acid as opposed to DNA. 

The main drawback of this work is that the fluctuation effects have been 
treated on a harmonic level which a priori assumes that they are small. 
However what saves us here is that the mean-field interactions decay ex- 
ponentially whereas the fluctuation interactions decay algebraically. This 
is different than in the case of a simple, i.e. unconnected, Coulomb fluid 
such as an inhomogeneous electrolyte. In that case both the mean-field as 
well as the fluctuation interactions decay exponentially but the mean-field 
interactions have a range twice as long as the fluctuation interactions |11]. 
The case of polyelectrolytes is in this respect fundamentally different: the 
range of fluctuation interactions is larger then the range of mean-field inter- 
actions and the use of harmonic approach to fluctuations is probably better 
grounded in this case. 
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